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SHORT METHOD OF ELLIPTIC FUNCTIONS. 



BY LEVI W. MEECH, A. M., HARTFORD, CONN. 

{Continued from page 136 ) 

When the last modulus is so small that it may be taken as 0, the denom- 
inator, |/(1 — e,? sin'tf,,), evidently reduces to its first term 1 ; whence F(0,0 n ) 
= d n ; where d n denotes the corresponding limit-amplitude. Then, 

F{e, 6) = n !(l+e»H(l+e 00 ) i(l+e w ) . . . f 
Let 

A = (1+0(1+0(1^)... = ^ . ^ . T ^... = ^-f^. 

Also let n denote the number of approximations, or of factors in the nume- 
rator; then 

F(e, d) = A.(6 n -+-2 n ) = Ax. 
When 6 = Jjt, the formula of amplitude gives 6' — n; 0" — In; 6'" 
= 4w; etc. Hence, for the quadrantal value, 

F(e, Jff) = A§;r. 
To aid in judging of the approximation, we take e=sin22^-°, or 6=cos22J° 
then c° = sin 2°16'03" = .03956, 6° = cos 2°16'03" = .99922, 

e oo = s i n o°13'28" = .00039, b°° = cos0°13'28" = .99999991, 
e«*> = .000 000 038 35+ b m = 999 999 999 999 999 26. 

Again, let e = sin 45°, or e 2 = \ — 0.5; then in common logarithms, 



Log e, sin 45° = 9.84948500 

" e°, sin 9°52'45" = 9.23444863 

" e 00 , sin 0°25'40' =7.87330122 

" e 000 ,sin 0°0'2" 877= 5.14455457 

£(e«*>) 2 = .000 000 000 048 65 



Log b, cos 45° = 9.84948500 

" 6° sin80°07'15" =9.99351181 
" b 00 sin 89°34'20" = 9.99998788 
" 6 000 sii, 89°59'57".1= 9.99999999 
" A == 0.072 007 344 8. 



r o = 2e oo _ (l=^)sin 2r+J(L^) 2 sin4^- 
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This last series is derived from the formula of amplitudes by the well 
known trigonometric development. When is less than 45°, the last writ- 
ten term can evidently be omitted, and the result will still be exact to ten 
decimal places ; that is (1 - 6 00 )-f- (1 + b w ) = Je 00 2 - itan%>°. And by re- 
duction, A = (- r |* . T = =, - ; = r^, 

' U + 6/ l-f-i/6 v /[(cos^+l). 7 /(cosl^)]' 

F{e, 0) = A.l(6 m — Je 002 sin2fl 00 ) = ^K^ 00 — |tan%>°sin 40°). 
When e 2 is greater than \, we may use the formulas of ascending moduli, 
Sect. IV, and the third approximate value or ff" will still be exact to ten 
decimals. 

Note 1. Another new scale of amplitudes will be developed in Section 
VII. 

Note 2. By shorter process from Sect. VIII, let W— 0-\- 2w+w; tanu 
= j/^.tan ; tan v — b.t&n ; then 

dO _J_ 1 dW 

,/(l+e 2 sin 2 0) "" 1 + b •' 1 + b° ' |/(1 — e'^sin 2 ?'')' 
Note 3. Again, as derived from Jacobi, let the common logarithm of q, 
that is, Log q = — 0.434294481 9 X nA'^A; where A' ■+- A = F(b, %iz) 
-*-F(e, £tt); also Log 0.434294482* = 0.1349341840; and by Sect. VII; I 

a ,r l—q—q 4 +2q 9 —q™— 4 W—0 , , n e 3 

J= Vo.- — a — 34 ! * — 3^ ; tan — -— == J.tanP: e —- ^- ; 

1+5— g 4 — 2g 9 — 5 16 + ' 2 ' {1 + Jy 

dd 1 d?P" 



then 



l/(l— e 2 sin 2 #) 1+22" p/fl— e 02 sin 2 r)' 



V. Integration by Series. As shown by analytical trigonometry, 
if e = sin <p , then <IF or 

d£_ dO <W 

A ~~ i/(l— f?+ie 2 cosW) " _J co7 2 |^i/( 1 H-tah 2 !^" V-l)(i +tanfy>rW = I)- 
Expanding and integrating 

-- 2 «[M+ra te "'I+il t »1+---] 

-T !! I>1+ -]+* 

In like manner, for the function -B, of the second species, 

dE — cos 2 | ? > ) /[(l+tan 2 ^.£2V- 1 )(l +tan , ^.ff-2V-l)]d0. 
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Expanding and integrating as before, 

E = eos 2 ^ { ^[l+(J) 2 tan 4 | f +(|) 2 tan 8 ip+(^) 2 tan 12 ^+ . . . .] 
+ sin 2d\_\tstf\y +^tan%+ T ^tan 10 ip + . . . .] 
— Jsin 4(9[|tan 4 if>+^tan 8 ^+ . . . .] 
+|sin 6d[ r \tan 6 i<p+ . . . .]+ete. } . 
When d = \tc the preceding series for F and E reduce to their first terms. 
Having computed these, the remaining terms can be accurately derived from 
each other by the following process from Legendre, (Fonotions Elliptiques,) 
Vol. 1, p. 276. 

F = P j- d6 2 . , m = Ad— A j sin 2d + A „ sin40— A 3 sin60 + . . . 
J V(l — e 2 sm 2 #) x * i < 

Differentiating once, and again, 

dF 1 



dd i/(l— e 2 sin 2 0) 
e 2 sin 2d 



A— 2A 1 cos 2d+AA 2 cos Ad— 6^ 3 cos 6 (?+ . 
-= 4J.jsin 20— 16A 2 sin40+36J. 3 sin60— . 



2/(1 — eWi 

Multiplying the last equation by 2(1 — e 2 sin 2 #) -t- e 2 , or by (2-^e 2 ) — 1-f 
cos 2d, and the equation preceding by sin2#, to render the left hand members 
identical, we equate the coefficients of sin 2d, sin Ad, .... in the right hand 
members; thus, we shall find, 

2X3^l 2 = 4JL 1 [(2-4-e 2 )— 1]— A, 
3X5^3 = 16J. 2 [(2-f-e 2 )— 1]— 1X3J. X , 
4X74 4 = 36^ 3 [(2^e 2 )-l]— 2X5^ 2 , 



5X9.4 5 == 64^ 4 [(2^e 2 )—l]— 3X7.A 



3' 



Again assuming for E and differentiating, we have, 

E =fdd i /(l—e 2 sin 2 d)=Bd+B 1 sm2d—B 2 sin 40+.B 3 sin60— . . . 

^= y(l— e 2 sin 2 0) = B+2B 1 eos 2d—AB 2 cosAd+QB 3 cos6d—... 

\ehm2d _ 4£ . 2 0_ 16 r sin 40+36.B„sin 60— . . . 
V(l— e*sin*0) 1 2 i s 

Multiplying the preceding derivative of F by ^e 2 sin 2d, to make the left 
hand member identical with this last equation ; than equating the coeffi- 
cients of sin 2d, sin Ad, .... in the right hand members, 

4*! = \<?{A-2A % \ 
16_B 2 = ie\A 1 — SA 3 ), 
36i? 3 = ie 2 (2^2-4A),.... 
An*B n = ^\Xn-\)A n _ r -{n+l)A n+1 \. 
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The equation for A l is easily found by substituting the series for F and 
for E in the identical equation, 

dE = Add = A*d6+ A = (1— V+JAos20)dF. 

Whence B = (l-^jA-^A 1} or A, = a(|_i)-|b. 

Had we developed directly by the binomial theorem, and then integrated, 
we should have found, for quadrants, 

^w-^-S[i + @V + ^S)W(j^)V + ...]. 

Here A = — e 2 — | — ). 

cteV e / 

The exponential developments have just given, with c = sin <p, 
A = c ^[ 1 +i tan4 ^+(t) 2tan8 ^+(A) 2 ^n 12 |?'+ • • •], 

£= cos%> [l+|tanH?+a) 2 tanH?'+(T 1 6) 2 tan 12 ^+---]- 
Lastly, to develope after a single decrease of modulus, let t be determined 
by the relation j/6 = ]/(cos ^>) = (1 — t)-t-(l-±-t), or 

e = siny>; e° = sin (p° = tan 2 -^; t = tan Jf° = Ve n °. 
For reference, F(e, ^n) = (1-f e°).F(e , f;r), l+e° == (l+Q»-g-(l+^) 

= (l+i) 2 cos 2 ^ . 

A = (l + ^.[l + (i)¥ + (f)¥+(^) 2 P + ( T 3 ^)^ 6 +(^) 2 ^+. • .]• 
When e = .9851714 = sin 80°07'15", f = 45°, t = tan22£°; ^ftfe 
x t K = .000 000 044. The first omitted term is .000 000 000 0318 ; so that 
the written terms would give the value of A exact to ten places. 

The common logarithms of the series of squares of binomial coefficients 
are, Log (J) 2 =1.39794 00087; Log (f ) 2 =1.14806 25354; 
Log tftf = 2.98970 00434 ; Log ( T 8 ^) 2 = 2.87371 61496 ; 
Log {J^f = 2.78220 11684 ; Log <J^f = 2.70662 40466. 

VI. Inverse Method oe Integration. Eeferring to the close of 
Section II, let the longest side of the spherical triangle o = 90°. The 
other two sides a, b, or 0, 0', termed complementary, are then connected by 
the equations, 

F{e, d)+F[e, 6') = F(e, %n); or A.{x+x') = A.\n; 
= cos cos 0'— sin sin #V(1— e 2 ). 
If J = \/(l— e 2 sin 2 0), e sin = 4/(1— J 2 ), ecos d = i/{A 2 — 1+e 2 ); and 
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so for 0', A'. Multiply the preceding equation by e 2 , substituting in terms 
of A, A', and reducing, since A' = \/{\ — e 2 sin 2 #') 

A A' = t/(1— e 2 ) = 6. 
This constancy of the complementary product leads to the preliminary 
assumptions, A — y'J./cos 20, 4' = ^6. /cos 20'. Here / must denote 
such a function of cos 2d that the product /cos 20./cos 2d' = 1 ; or one 
function is equal to the reciprocal of the other. And since 2dF — Ad(2x) 
= d(2d) ■+■ J, and by the Calculus, cZ(2a;) is — e/ cos 2a; -f- 1/(1 — cos 2 2a;), it 
follows that /cos 20 is /'cos 2a;. Changing then to 2a;, since 2a;' = ic — 2a;, 
as above indicated, we see that /' must denote such a function that /'cos 2a; 
is the reciprocal of /'cos (iz — 2a;). This condition is fulfilled by assuming 

■j/(l — e 2 sin 2 0) l+^iCos 2a;+a 2 cos 4a;+a 3 cos 6a!+a 4 cos Sx-\-. . . 

y 6 ~~ 1 — t^cos 2a;-|-a 2 cos ^ x — a z cos 6a;+a 4 cos 8a; — . . . 

Or adding and subtracting the identical equation 1 = 1, and dividing 
the one result by the other; also denoting the quotient by 8, 

o A — i/b a t cos 2a;-f a 3 cos 6a;-fa 5 cos 10a;-|- • • • 

A-\-y/b l-|-a 2 cos4a;-fa 4 cos8a;+ . . . 
Multiplying both sides by the right hand denominator, differentiating 
and dividing by — dx; since dd = AAdx, 

(l-\-a 2 GOS Ax+a^s 8x-\-. • •) = 2a x sin 2a; -t-6a 3 sin 6a; +. .. 
In the last two equations, omitting all terms of the series after a 3 , let us 
substitute the correlative values, firstly x = 0, d =z 0, or A =z 1. The for- 
mer equation gives, since, by Section V, -j/6 = (1 — t)-^(l-\-t), 
c&i — fa 2 + a 3 — ta i = t. 
But the derivative or latter equation vanishes. Secondly, the correlative 
values x = 45° or 2a; = 90°, A = A' — \/b, cause the former equation 
to vanish, while the latter becomes, since e sin d = -j/(l — A\ etc., 
a 1 +iA(l—b).a^—Sa s —iA(l~b)a i = +iA(l—b). 
Approximately, omitting a &) a 4 , and taking A, b in terms of t from Sec- 
tion VI, we find from the last two equations, 

a l =t + tf + ...,a a = &+... 
Again differentiating, and substituting in place of x and d or A, as before, 
we find a third equation; from which and the two former, omitting a 4 , we 
find by elimination, exact to another power of t, 

a x = t + -gt -\- ^-g^t + . . ., a 2 = -gt -f- j-gt + • • •, a 3 z=: - 3TS^ + • • • 

By repeated differentiations, other equations can be found in the same 
way, and the approximation carried to any extent, disclosing this remark- 
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able law; let t = (1 _,/&)-(l + ,/&); q = \t + ft? + J&P + tfftp* 

+ tHIt2 <17 + • • -J *«» «i = H o a = 2 9 4 > « 8 = V, «4 = 2 S 16 > «s 
= 2<f 5 , a 6 — 2g 36 , .... This has been fully demonstrated in a different 

way by its illustrious discoverer, Jacobi. He has given in the twenty-sixth 

volume of Crelle's Journal an extended table of the values of q for every 

6' of f . In this, and in other ways may be demonstrated one of the most 

important results, the inverse relation, 

T/(l-e a sin 8 g) __ 1 +2q cos 2a;+2g 4 cos 4a?+29 9 cos 6x+ ... 
Vb " 1 — 2q cos 2x+2q*cos 4% — 2g 9 cos 6a; + . . . 

In the next place, for determining the integral x through cos 2%, by 
reversion of series, we take the second form, found before, 

a J — f/b 2 gcos 2x 4- 2<fcos Qx -|- . . . . 

" " J+\/b ~ 1 +2g 4 cos 4«-r2g 19 cos 8a:+.77~ 

Multiplying both sides by the right hand denominator, and changing 
cos 4a;, cos 6a;, etc., to powers of cos 2x, also omitting the sixteenth and 
higher powers of q. 

/S(l— 2g 4 +4g 4 cos 2 2a;) = (2q—6q 9 )cos2x(l+4q s cos i 2x). 
Dividing both members by the right hand factor, omitting q n , 

S[l— 2q i +4{q*— g^cos^a;] = (2q—6q 9 )cos 2x. 
This equation, being a common quadratic, is easily resolved through an 
auxiliary arc v; the well known form of the result will be, when P and Q 
are functions of q only, sin v = PS • cos 2a; = Q tan ^y. 

To determine P and § independently of q, we first elimiminate v, since 
sin«z=; 2tan^v-r- (l+tan 2 ^»), and in the first and third derivatives of the 
result, we make x = 45°, and J — \/b, which gives, if sin <p — e, cos <p 
— b; t =z (1 — •j/6)-=-(l+>' b), the practical solution, 

f) . I r ^§ I 2 9_/2 16 5/6 

The resulting integral a; should be exact to <? 12 , that is to ten decimal 
places when <p is less than 70°, or exact to seven decimals when <p is less 
than 83°. The auxiliary S varies from the limit — t to -\-t ; so that v may 
be positive, 0, or negative ; and 2a; greater or less than 90° or \x. When 
(p is near to 90°, the value of x may be found from the standard equation 
by the process of trial and error, to any degree of accuracy. Finally 
F — A x x, 

(To be concluded in No. 1, Vol, V.) 



